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We investigate the GRB formation rate based on the relation between the 
. spectral peak energy (Ep) and the isotropic luminosity. The E'p-luminosity re- 

! lation covers the energy range of 50 - 2000 keV and the luminosity range of 

^ ! 10^°-10^^ erg s~^, respectively. We find that the relation is considerably tighter 

^ ! than similar relations suggested previously. Using i?p-luminosity relation, we es- 

'. timate the luminosity and the redshift of 684 GRBs with the unknown distances 



O 

:^ 

6 



and derive the GRB formation rate as a function of the redshift. For < 2; < 2 
, the GRB formation rate is well correlated with the star formation rate while 
it increases monotonously from 2; ~ 2 out to z ~ 12. This behavior is consis- 
tent with the results of previous works using the lag-luminosity relation or the 



j§ ! variability-luminosity relation. 



^ ! Subject headings: gamma rays: bursts — cosmology: early universe 



1. Introduction 

Many ground based telescopes have tried to detect optical afterglows of GRBs and 
to measure their redshifts by using the spectral absorption and/or emission lines of the 
interstellar matter in the host galaxy. However, the number of GRBs with measured redshift 
is only a fraction of all GRBs detected with BATSE, BeppoSAX, HETE-II and INTEGRAL 
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satellites; we have still only about 30 GRBs with the known redshifts. The most of them 
occur at the cosmological distance, and the current record holder is GRB 000131 at z = 4.5 
(Andersen et al. 2000). According to the brightness distribution of GRBs with the known 
redshifts, the above satellites should have already detected much more distant GRBs, such as 
at z ~ 20 (Band 2003). If we can establish a method for estimating the intrinsic brightness 
from the characteristics of the prompt gamma-ray emission, we can use the brightness of 
the GRB as a hghthouse to determine the unknown redshifts of majority of GRBs, which 
enables us to explore the early universe out to z ~ 20. 

Using the geometrical corrections of the collimatcd jets. Frail et al. (2001) and Bloom et 
al. (2003) revealed that the bolometric energies released in the prompt emission are tightly 
clustering around the standard energy of ~ 1 x 10^^ erg. Thus, the explosion energy of GRBs 
can be used as a standard candle as the supernovae. However, the apparent brightness of 
GRBs strongly depends on the jet opening angle and the viewing angle. To use the GRB as 
the standard candle, we need to correct such effects. 

Several authors tried to establish a method for estimating the isotropic luminosity from 
the observed GRB properties. Using the lightcurves of the prompt gamma-ray emissions, 
some pioneering works have been done; the variability-luminosity relation reported by Fen- 
imore & Ramirez- Ruiz (2000), which indicates that the variable GRBs are much brighter 
than the smooth ones. The spectral time-lag, which is the interval of the peak arrival times 
between two different energy bands, also correlates with the isotropic luminosity (Norris et 
al. 2000). These properties based on the time-series data might be due to the effect of the 
viewing angle from the GRB jet (e.g., loka & Nakamura 2001; Norris 2002; Murakami et al. 
2003). 

On the other hand, based on the spectral analyses with the K-correction (Bloom et 
al. 2001), Amati et al. (2002) found the correlation between the isotropic-equivalcnt energy 
radiated in GRBs and the peak energies Ep, which is the energy at the peak of z/iv spectrum. 
Atteia (2003) suggested a possibility of the empirical redshift indicator, which is based on 
the Ep and the arrival number of photons. 

Applying these luminosity indicators to the GRBs with the unknown redshifts, their 
redshifts can be estimated from the apparent gamma-ray brightness. As a natural application 
of the obtained redshift distribution, the GRB formation rates are discussed by several 
authors (e.g., Fenimore & Ramirez-Ruiz 2000; Norris et al. 2000; Schaefer et al. 2001; Lloyd- 
Ronning et al. 2002; Murakami et al. 2003). Especially, using mathematically rigid method 
(Efron & Petrosian 1992; Petrosian 1993; Maloney & Petrosian 1999), Lloyd- Ronning et al. 
(2002) have estimated the GRB formation rates from the variability-luminosity relation. 
These works give basically the same results; the GRB formation rate rapidly increases with 
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the redshift at < 2; < 2, and it keeps on rising up to higher redshift {z ~ 12). The GRB 
formation rate did not decrease with z in contrast with the star formation rates (SFRs) 
measured in UV, optical and infrared band (e.g., Madau et aL 1996; Barger et aL 2000; 
Stanway et aL 2003). 

In this paper, we estabhsh a new cahbration formula of the redshift based on the Ep- 
luminosity relation of the brightest peak of the prompt gamma-ray emission from 9 GRBs 
with the known redshifts. Importantly, the uncertainty of our formula is much less than 
those of the previous works (lag, variability) as shown in section 3. Applying the obtained 
calibration, in section 4, we estimate the redshifts of 684 GRBs with the unknown redshifts. 
We then demonstrate the GRB formation rate out to 2; ~ 12 and the luminosity evolution 
using the non-parametric method (Efron & Petrosian 1992; Petrosian 1993; Maloney & 
Petrosian 1999; Lloyd- Ronning et al. 2002). We emphasize that the GRB formation rate 
derived by us is based on the spectral analysis for the first time, and its uncertainty is well 
controlled in a small level. Throughout the paper, we assume the flat-isotropic universe with 
flm = 0.32, Oa = 0.68 and Ho = 72 km s'^Mpc"^ (Bennett et al. 2003; Spergel et al. 2003). 



We used 9 GRBs (970508, 970828, 971214, 980703, 990123, 990506, 990510, 991216, 
and 000131) in the BATSE archive with the known redshifts, and focused on the brightest 
peak in each GRB. We performed the spectral analysis with the standard data reduction for 
each GRB. We subtracted the background spectrum, which was derived from the average 
spectrum before and after the GRB in the same data set. 

We adopted the spectral model of smoothly broken power-law (Band et al. 1993). The 
model function is described below. 



for E < {a—P)EQ and E > {a—P)EQ, respectively. Here, N{E) is in units of photons cm~^s"^keV 
and Eq is the energy at the spectral break, a and /3 is the low- and high-energy power-law 
index, respectively. For the case of /3 < — 2 and a > — 2, the peak energy can be derived as 
Ep — {2 + a)Eo, which corresponds to the energy at the maximum flux in i/Fj, spectra. The 
isotropic luminosity can be calculated with the observed flux as L = AirdlF^, where di, and 
F^ are the luminosity distance and the observed energy flux, respectively. 



2. Data Analysis 




(1) 
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3. Luminosity relation 

In figure 1, we show the observed isotropic luminosity in units of 10^^ erg s~^ as a 
function of the peak energy, Ep{l + z), in the rest frame of each GRB. For one of our sample 
(GRB 980703), a lower limit of Ep{l + z) is set because the spectral index /3 > —2. The 
results of BeppoSAX reported by Amati et al. (2002) are also plotted in the same figure. 
Here, we converted their peak fluxes into the same energy range of 30 - 10000 keV in our 
analysis with their spectral parameters. There is a good positive correlation between the 
Ep{l + z) and the L. The linear correlation coefficient including the weighting factors is 0.957 
for 15 degree of freedom (17 samples with firm redshift estimates) for the log[£^p(l + z)] and 
the log[L]. The chance probability is 1.85 x 10~^. When we adopt the power-law model as 
the £^p-luminosity relation, the best-fit function is 



EJl + z 



1.94±0.19 



— = (4.29 ± 0.15) X 10-^ ^\ ' (2) 

1052 ergs ^ ^ L 1 keV J ^ ' 

where the uncertainties are 1 a error. This relation agrees well with the standard synchrotron 
model, L(xE^ (e.g., Zhang & Meszaros 2002; Lloyd et al. 2000). 



4. Redshift Estimation and GRB Formation Rate 

The £'p-luminosity relation derived in the previous section seems to be much better 
estimator of the isotropic luminosity compared with the spectral time-lag and the variability 
of GRBs (Norris et al. 2000; Fenimore & Ramirez-Ruiz 2000; Schaefer et al. 2001) since the 
chance probability is extremely low. In this section, using the £^p-luminosity relation, we 
try to estimate the isotropic luminosities and the redshifts of the BATSE GRBs with the 
unknown redshifts. 

We first selected about 1000 GRBs in a class of the long duration of Tgo > 2 sec 
detected by BATSE, and extracted the brightest pulse in each GRB. We performed the 
spectral analysis for these peaks with the same method described in section 2. Once we 
obtained observed energy-flux and Ep at the observer's rest frame, we can estimate the 
redshift with equation (2) and the luminosity distance as a function of the redshift. We 
could not determine Ep for ~ 5 % of samples because oi (3 > —2 within 1 a uncertainty, 
and excluded them. After setting the flux limit of Fumit = 1 x 10^^ erg cm~^s~^, which is 
based on the dimmest one (GRB 970508), on the data set of {z,L) plane, 684 samples are 
remained. In figure 2, we show the sample distribution in [z, L) plane with the truncation 
by the flux limit. 

The normalized cumulative luminosity function at each redshift bin is also shown in 
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figure 2 . The shape of these luminosity functions are similar to each other, but the break- 
luminosity seems to increase toward the higher rcdshift. This indicates that a luminosity 
evolution is hidden in the {z, L) plane in figure 2. Therefore, we have to remove the effect 
of the luminosity evolution from the data set before discussing the real luminosity function 
and the GRB formation rate. This is because the univariate distribution of the redshift and 
the luminosity can be estimated only when they are independent of each other. 

We used the non-parametric method (Lyndcn-Bell 1971; Efron & Petrosian 1992; Pet- 
rosian 1993; Maloney & Petrosian 1999), which were used for the Quasar samples and first 
applied to the GRB samples by Lloyd- Ronning ct al. (2002). The total luminosity function 
^{L,z) can be written ^{L,z) = p{z)(l){L/gk{z),as)/gk{z) without the loss of generality. 
Here, each function means the luminosity evolution gk{z), the density evolution p{z) and the 
local luminosity function (f){L / g}.{z) , cts), respectively. Although the parameter represents 
the shape of the luminosity function, we will ignored the effect of this parameter because 
the shape of the luminosity function is approximately same as shown in figure 2. We do not 
mention more about the method, but the details are found in Maloney & Petrosian (1999). 

Assuming the functional form of the luminosity evolution as gk{z) = (1 -|- z^ , which 
is also used in Maloney &: Petrosian (1999) and Lloyd-Ronning et al. (2002), we convert 

the data set of {z.L) to {z,L'), where L' = L/gk{z). We search the best value of k giving 
the independent set of (z, L') within the significance of 1 o" error, and it is found to be 
9k{z) = (1 + zY'^^^^^''^^ . Once we obtain the independent data set of the {z,L'), we can 
generate the cumulative luminosity function ip{L') with a simple formula (see equation (14) 
of Maloney & Petrosian 1999). In figure 3, we show the cumulative luminosity function 
calculated by the formula. This function is approximately described by '4>{L') cx L'~°-^ for 
0.1 < L'52 < 1 and V'(^') oc L'-^ ^ for 1 < L'52 < 50, where L'52 = L'/IO^^ ergs. 

We can also obtain the cumulative number distribution ip{z) as a function of z. The 
differential form of the function is useful for the purpose of comparison with the star forma- 
tion rates in other wave bands. We convert '(p{z) into the differential form with the following 
equation. 



where the additional factor of (1+^) comes from the cosmological time dilation, and dV{z) / dz 
is a differential comoving volume. In figure 4, we show the relative comoving GRB rate p{z) 
in unit proper volume. 




(3) 
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5. Discussion 

We have investigated the spectral property of the brightest peak of each GRB with the 
known redshifts, and have found a fine correlation between the peak energy Ep{l + z) and the 
isotropic luminosity. While this correlation against a smaller sample has been pointed out 
by some authors (e.g., Amati et al. 2002; Atteia 2003; Schaefer 2003a,b), we have succeeded 
in combining the results of BeppoSAX and BATSE to describe equation (2). 

Using the new £^p-luminosity relation, we have estimated the redshifts of 684 GRBs 
with the unknown redshifts. We found the existence of the luminosity evolution gk{z) = 
(1 + 2;)^-^^='=° °^ for GRB samples as shown in figure 2. The null- hypothesis of the luminosity 
evolution is rejected about 9a significance. Lloyd-Ronning et al. (2002) also suggest the 
presence of the luminosity evolution as gk{z) — {l + zY''^^^-^. These two values are consistent 
with each other. The luminosity evolutions are found in other objects. For example, Caditz 
& Petrosian (1990) and Maloney & Petrosian (1999) estimated the luminosity evolution of 
the QSO samples as gk{z) = (1 + and (1 + z)^-^^^. respectively. Based on the observation 
of Subaru Deep Field and a photometric redshift estimation for K'-hand selected galaxy 
samples, Kashikawa et al. (2003) found the strong luminosity evolution in the rest UV band 
of their galaxies. 

The form of the cumulative luminosity function is independent of z except for the 
break luminosity, which changes with z. We propose that figure 3 might include important 
information on the jets responsible for the prompt 7-ray emissions and a distribution of 
their opening angles. Consider an extremely simple model for a uniform jet with an opening 
half-angle 9j and a constant geometrically-corrected luminosity Lq, which is viewed from an 
angle of 9y. Then, in a crude approximation the luminosity L is given by 

^ r 2Loef for 9, < 9, 

\ 2L,{9y9l) for 9, > 9,. ^ > 

For the case of 9,^ > 9j, L is proportional to d^, where 6 = [7(1 — f3 cos9,ij)]^^ oc so 
that the luminosity has the dependence of 9^^ (loka & Nakamura 2001; Yamazaki et al. 
2002, 2003). The dependence of 9j is determined in order that two functions in equation (4) 
are continuously connected at 9^, = 9j. We also consider the distribution of 9j in the form 
f{9j)d9j oc 9~'^d9j when ^^min < 9j < ^^max- Then, easy calculations show that in the case of 
q < 5/2, we have 

« I i(,-3,/. fo. < i < 2Lo9- . ■ 

This is a broken power low with the break luminosity 2Lq9^^. Then if ^^L-^o oc gk{z) with 



- 7- 



q — 0.6, we can roughly reproduce figure 3. This suggests that either the maximum opening 
half angle of the jet decreases or Lq increases as a function of the redshift. 

This work is the first study to generate the GRB formation rate with £'p-luminosity 
relation. The result indicates that the GRB formation rate always increases toward z ~ 12. 
This is consistent with the previous works using the GRB variability (Fenimore & Ramirez- 
Ruiz 2000; Lloyd-Ronning et al. 2002) and the spectral time-lag (Norris et al. 2000; Schaefer 
et al. 2001; Murakami et al. 2003). Quantitatively, GRB formation rate is proportional to 
(1 + z)^^^ foT z < 1 and to {1 -\- z^ for z > 1. On the other hand, the SFR measured in UV, 
optical, and infrared is proportional to (1 + z)'^^^ for z < 1 (e.g., Lilly et al. 1996; Cowie 
et al. 1999; Glazebrook et al. 2003) and to (1 + z)-^-^ ioi z > 1 (e.g., Madau et al. 1996; 
Barger et al. 2000; Stanway et al. 2003; Kashikawa et al. 2003). Therefore, we find that the 
ratio of the GRB formation rate to the SFR evolves along redshift, say oc (1 + z)^~^. 

Recently, it has been strongly suggested that the long duration GRBs arise from the 
collapse of a massive star (Galama et al. 1998; Hjorth et al. 2003; Price et al. 2000; Uemura 
et al. 2003; Stanek et al. 2003). Hence our result imphes that either the formation rate of the 
massive star or the fraction of GRB progenitor in massive stars at the high redshift should 
be significantly greater than the present value. However, if the SFR rapidly increases along 
redshift as suggested by Lanzetta et al. (2003), the fraction of GRB progenitors does not 
change so much. 

The existence of the luminosity evolution of GRBs, i.e., gk{z) — (1 -|-2;)^-^^, may suggest 
the evolution of GRB progenitor itself (e.g., mass) or the jet evolution. Although the jet 
opening angle evolution was suggested (Lloyd-Ronning et al. 2002), in our extremely simple 

model, either the maximum jet opening angle decreases or the jet total energy increases. In 
the former case the GRB formation rate shown in figure 4 may be an underestimate since 
the chance probability to observe the high redshift GRB will decrease. If so, the evolution 
of the ratio of the GRB formation rate to the SFR becomes more rapid. On the other hand, 
in the latter case, GRB formation rate of figure 4 gives a reasonable estimate. 
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analysis tools. This research is supported by Grants-in-Aid from the Japanese Ministry of 
Education, Culture, Sports, Science and Technology (15740149, 14204024, 14047212 and 
12640302). R.Y., A.K.I, and K.I. are supported by JSPS Research Fellowship for Young 
Scientists. 
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Fig. 1. — The £'p-luminosity relation. The open squares are our new results with BATSE. 
The results of BeppoSAX (Amati et al. 2002), which are converted into the energy range of 
30 - 10000 keV, are also shown as the filled squares and the cross points. The sohd line is 
the best-fit power-law model for the data. The linear correlation coefficient is 0.96 for 15 
degree of freedom. 
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Redshift 



Fig. 2. — The distribution of luminosity vs. redshift derived from the E^p-luminosity relation. 
The truncation of the lower end of the luminosity is caused by the flux limit of Fumit = 
1 X 10~^ erg cm^^s""^. The inserted figure is the cumulative luminosity function in the 
several redshift ranges. The luminosity evolution exists because the break-luminosity increase 
toward the higher redshift. 
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Fig. 3. — The cumulative luminosity function of L' — Lj g^iz). This function is equivalent 
to the present distribution at z = because the effect of the luminosity evolution g^iz) is 
removed. 



-13- 



T — r 




L I I I I I I I I l_ 

12 5 10 



1+z 

Fig. 4. — The relative GRB formation rate normalized at the first point. The solid line is the 
result based on the best fit of £^p-luminosity relation and two dotted lines indicate the upper 
and lower bounds caused by the uncertainty of ii^p-luminosity relation. These dotted lines 
are also normalized and superposed on the best result at < 2; < 1 with the least-square 
method. The error bars accompanying open squares represent the statistical uncertainty of 
each point. 



